home *** CD-ROM | disk | FTP | other *** search
/ NASA Climatology Interdisciplinary Data Collection / NASA Climatology Interdisciplinary Data Collection - Disc 4.iso / software / grads / lib / functions / colorbar.gf < prev    next >
Encoding:
Text File  |  1994-02-17  |  1.7 KB  |  96 lines
  1.  
  2. *
  3. *  Script to plot a colorbar
  4. *
  5. function colorbar (shdinfo)
  6. *
  7. *  Check shading information
  8. *
  9.   if (subwrd(shdinfo,1)='None') 
  10.     say 'Cannot plot color bar: No shading information'
  11.     return
  12.   endif
  13. * 
  14. *  Get plot size info
  15. *
  16.   'query gxinfo'
  17.   say result
  18.   rec2 = sublin(result,2)
  19.   rec3 = sublin(result,3)
  20.   rec4 = sublin(result,4)
  21.   xsiz = subwrd(rec2,4)
  22.   ysiz = subwrd(rec2,6)
  23.   ylo = subwrd(rec4,4)
  24.   xhi = subwrd(rec3,6)
  25.   xd = xsiz - xhi
  26.   say ' xs   'xsiz' 'ysiz' 'ylo' 'xhi' 'xd 
  27. *
  28. *  Decide if horizontal or vertical color bar
  29. *  and set up constants.
  30. *
  31. *  if (ylo<0.6 & xd<1.0) 
  32. *    say "Not enough room in plot for a colorbar"
  33. *    return
  34. *  endif
  35.   cnum = subwrd(shdinfo,5)
  36.  
  37.   ylo=0.5
  38.   xd=1.0
  39.   xhi=7.5
  40.   ysiz=11
  41.   xsiz=8
  42.   vert=1
  43.  
  44.   if (vert=1)
  45.     xl = xhi + xd/2 - 0.4
  46.     xr = xl + 0.2
  47.     xwid = 0.2
  48.     ywid = 0.5
  49.     if (ywid*cnum > ysiz*0.8) 
  50.       ywid = ysiz*0.8/cnum
  51.     endif
  52.     ymid = ysiz/2
  53.     yb = ymid - ywid*cnum/2
  54.     'set string 1 l 5'
  55.     vert = 1
  56.   else
  57.     ymid = ylo/2
  58.     yt = ymid + 0.2
  59.     yb = ymid
  60.     xmid = xsiz/2
  61.     xwid = 0.8
  62.     if (xwid*cnum > xsiz*0.8)
  63.       xwid = xsiz*0.8/cnum
  64.     endif
  65.     xl = xmid - xwid*cnum/2
  66.     'set string 1 tc 5'
  67.     vert = 0
  68.   endif
  69. *
  70. *  Plot colorbar
  71. *
  72.   'set strsiz 0.12 0.13'
  73.   num = 0
  74.   while (num<cnum) 
  75.     rec = sublin(shdinfo,num+2)
  76.     col = subwrd(rec,1)
  77.     hi = subwrd(rec,3)
  78.     'set line 'col
  79.     if (vert) 
  80.       yt = yb + ywid
  81.     else 
  82.       xr = xl + xwid
  83.     endif
  84.     'draw recf 'xl' 'yb' 'xr' 'yt
  85.     if (num<cnum-1)
  86.       if (vert) 
  87.         'draw string '%(xr+0.05)%' 'yt' 'hi
  88.       else
  89.         'draw string 'xr' '%(yb-0.05)%' 'hi
  90.       endif
  91.     endif
  92.     num = num + 1
  93.     if (vert); yb = yt;
  94.     else; xl = xr; endif;
  95.   endwhile
  96.